Quantum statistics of the collective excitations of an atomic ensemble inside a cavity 
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We study the quantum statistical properties of the collective excitations of an atomic ensemble inside a high- 
finesse cavity. In the large-detuning regime, it is found that the virtual photon exchange can induce a long-range 
interaction between atoms, which results in correlated excitations. In particular, the atomic blockade phe- 
nomenon occurs when the induced long-range interaction effectively suppresses the double atomic excitation, 
when the average photon number takes certain values, which makes the two nearest energy levels degenerate. We 
also show that quantum phase transitions occur in the indirectly-interacting atomic ensemble when the average 
photon number reaches several critical points. In this sense, the quantum statistical properties of the collective 
excitations are very sensitive to the change of the average photon number. Our model exhibits quantum phase 
transitions similar to the ones in the Lipkin-Meshkov-Glick model. Our proposal could be implemented in a 
variety of systems including cavity quantum electrodynamics (QED), Bose-Einstein condensates, and circuit 
QED. 

PACS numbers: 42.50.Nn, 42.50.Dv, 73.43.Nq 
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I. INTRODUCTION 

In quantum optics, photon statistics reflect the essential 
properties of the electromagnetic field [1]. Importantly, cor- 
related photon counting by the second-order correlation func- 
tion can characterize the very quantum nature of light, such 
as bunching and antibunching effects [2], as well as the pho- 
ton blockade [3, 4], which is also referred to as optical state 
truncation [5]. The quantum statistical approach for photon 
counting [6] is also applicable to other massive and massless 
bosons [7]. The collective excitations of an atomic ensemble 
could be regarded as an operational quantum memory [8, 9] 
and the ensemble behaves as a boson in the large limit 
with low excitations [10]. Therefore, it is expected that the 
quantum statistical approach can also work well for atomic 
collective excitations. Moreover, the quantum correlations of 
these excitations can also be responsible for double-excitation 
effects, such as the Rydberg blockade, where double excita- 
tion is strongly suppressed by the dipole-dipole interaction 
between highly excited Rydberg atoms [11-13]. 

The atomic blockade is similar to the Coulomb blockade, a 
typical mesoscopic phenomenon where a single electron pre- 
vents an electric current from crossing some confined nanos- 
tructure [14—17]. Similar blockade effects have been predicted 
and also observed in quantum optical system for photons [3, 4] 
and cold atoms [11-13, 18]. Recently, phonon blockade has 
been studied [19]. The blockade effect, whereby a single par- 
ticle prevents the flow [3, 4, 14-16, 18] or excitation of many 
particles, provides a mechanism for the precise manipulation 
of quantum states of microscopic quantum objects at the level 
of a single particle. In this sense, it is essential for the imple- 
mentation of single-particle-based quantum devices. The pho- 
ton blockade effect may have applications in single-photon 
sources, needed for the physical implementation of quantum 
cryptography protocols [20] . 

In this paper we consider quantum correlation effects for 



an atomic system. One of the correlation effects studied is 
the Rydberg blockade effect. We consider a similar atomic 
blockade effect using an indirect-interaction coupling, which 
is induced by some confined photons in a cavity rather than by 
dipole-dipole interactions between atoms, as in the Rydberg 
blockade. Physical properties of atomic ensembles can also 
be quantified via spin squeezing [21]. 

Specifically, we study the case where an ensemble of two- 
level atoms are coupled to a cavity field with a large detuning 
frequency. The photons in the cavity can induce excitation 
hopping among atoms, which form a collective excited state 
described by the number of excited atoms. We will consider 
the case where the number of excited atoms is similar to the 
difference between the numbers of excited atoms and unex- 
cited atoms. Furthermore, the variation of half of this differ- 
ence equals the variation of the number of excited atoms. 

Similar to the generic Coulomb interaction for the Rydberg 
blockade [11-13], the induced interaction by cavity photons 
is also a long-range interaction and results in inhomogeneous 
energy-level spacings. More specifically, the structure of the 
energy levels depends on the average photon number We find 
that there will be two degenerate energy levels at an integral 
multiple of 1 /2 for the average photon number. If the average 
photon number slightly deviates from an odd multiple of 1/2, 
these two degenerate levels will become nearly degenerate but 
far away from other energy levels. Hence, it is difficult for the 
atomic ensemble to transit from the nearly-degenerate levels 
to other levels. This shows that the double excitation requires 
higher energy, which is oflF-resonant to two single excitations. 
Therefore, the atomic blockade effect could occur. If we fur- 
ther change the average photon number, the pair of nearly- 
degenerate energy levels shifts far away from each other, but 
one of them could end up closer to a neighboring energy level 
which was far away from this pair before changing the aver- 
age photon number Thus, the occurrence of atomic blockade 
can be controlled by the average photon number in the cavity. 
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(a) 

FIG. 1 : (Color online) (a) Schematic of a cavity field of frequency 
a> coupled to an atomic gas consisting of N two-level atoms with 
energy-level spacing a>A. A driving field of strength Clj and fre- 
quency aid is applied to the atoms; (b) The coupling diagram of one 
of the two-level atoms in the cavity. Here, A is the detuning between 
the atomic level spacing oja and the cavity field frequency ai, namely, 
A = cl)a - (Ji, and cjd is the frequency of the drive. 



Meanwhile, a quantum phase transition (QPT) [22-24] oc- 
curs when the average photon number is a half-integer, for 
negative detuning (the difference between the atomic energy- 
level spacing and the frequency of the cavity field). This 
is partially due to the energy-level crossing under the above 
conditions. The ground state changes drastically around the 
critical points characterized by the average photon number. 
This QPT behavior is similar to that of the Lipkin-Meshkov- 
Glick (LMG) model [25], which was studied in the quantum- 
information-process context in, for example, Ref. [26]. In this 
sense, we can regard our system as a modification of the LMG 
model. However, the critical points in our system are average- 
photon-number-dependent. This provides a controllable way 
to manipulate the system between different phases. 

To characterize various correlation phenomena of the 
atomic collective excitation, such as the atomic blockade and 
sensitivity of the QPT [27-31], we introduce a generalized 
second-order coherence function by replacing the annihilation 
(creation) a (a ' ) operator of photons in the usual second-order 
coherence function of photons with the lowering (raising) /_ 
(7+) operator of the collective atomic excitations. We prove 
that the antibunching effect occurs near odd multiples of 1 /2 
for the photon number, which implies that the double atomic 
excitation is suppressed. We also find significantly different 



behaviors on either side of the critical points. 

This paper is organized as follows. In Sec. II, we describe 
the system based on the Dicke model [29, 32]. The effective 
Hamiltonian is given in terms of the collective excitation of 
the atomic ensemble, and the ground state is analyzed for dif- 
ferent APNs. In Sec. Ill we then coherently drive the atomic 
ensemble and derive the effective Hamiltonian near two criti- 
cal points «^ =1/2 and «'j = y - 1 /2. In Sec. IV, we introduce 
the generalized second-order coherence function and calcu- 
late the statistical properties of the excitations of the atomic 
ensemble in the cases with and without dissipation. We dis- 
cuss the atomic blockade effect and sensitivity of the QPT to 
the photon number in Sees. V and VI, respectively. Finally, we 
present our conclusions in Sec. VII. The explicit form of the 
parameters used in Sees. IV and VI are given in the appendix. 



II. QUANTUM CRITICALITY OF AN ATOMIC 
ENSEMBLE STRONGLY COUPLED TO A CAVITY FIELD 

A. Model and Hamiltonian 

As shown in Fig. 1, the system we consider consists of an 
ensemble of atoms confined to a single-mode cavity of fre- 
quency a>. The cavity field is described by the annihilation 
(creation) operator a (a ' ). This model can be implemented in 
a variety of systems including cavity QED [33], Bose-Einstein 
condensates [34], and circuit QED [35]. 

Our model is described by the Dicke Hamiltonian [29, 32, 
35^0] (hereafter, we take Ji - I), 

H, = uja^a + ^ y o-['> + -^f + acrf) , (1) 

under the rotating-wave approximation. Here, we use the 
PauU matrices cr9' = \e)jj{e\ - \g)jj{g\, crV = \e)jj{8\, and 
o"*-' = \g}jj{e\ to describe the atomic transition of the jth atom 
with energy-level spacing oja, where \e}j and \g)j are the ex- 
cited and ground states of the jth atom, respectively. 

For an atomic gas with size smaller than the wave- 
length [29, 31], we assume that all the atoms are located 
near the origin and interact with the cavity field at the ho- 
mogeneous coupling rate go/ y/N. Here, the factor in 
the denominator of the coupling strength originates from the 
fact that the coupling strength is inversely proportional to the 
square root of the volume of the cavity field 1 / Vv. The vol- 
ume V is approximately equal to the total volume occupied 
by the atoms, which is times the volume of a single atom. 
Hence we can write the factor y/N explicitly in the coupling 
strength. 

We would like to point out that, the superradiant phase tran- 
sition based on the Dicke model in a real atomic system does 
not exist due to the inclusion of electromagnetic vector poten- 
tial A- term [29, 37, 41, 42]. However, the following argu- 
ments about QPT are based on the LMG model [25], which 
will be derived from the above Dicke model, even including 
the term. The similar A^ term {V^ term) in circuit QED 
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system will not influence the Hamiltonian significantly, ex- 
cept for just a little shift of the critical point [43]. 

The atoms we consider are largely detuned from the fre- 
quency Lo of the cavity field; namely, the detuning A{= <jJa - 
cS) is much larger than the corresponding coupling strength 
gdl '\/N, that is, |A| » |go/ "V^j- In this case, one can use the 
Frohlich-Nakajima transformation [44, 45] (or adiabatic elim- 
ination method), to obtain the effective Hamiltonian, 



+ \w J (cr^'V^i^' + ^r'i'V/^^), (2) 

where W - g^/{NA) is the strength of the effective interac- 
tion among the atoms, which is induced by the virtual photon 
exchanges. The form of the Hamiltonian is very similar to 
the dipole-dipole interaction of atoms in free space. We note 
that the Frohlich-Nakajima transformation is equivalent to the 
approach based on the adiabatical elimination and some per- 
turbation theories [46]. Furthermore, the photon number a' a 
becomes a conserved number 



B. Symmetric Hilbert Space and the LMG model 

We now describe the Hilbert space of the symmetric exci- 
tation. The Hilbert space of two-level atoms is spanned by 
2^ basis vectors {\gi), |e/)} with / = 1 , 2, ■ ■ ■ , A^. In the present 
case, all the atoms have identical transition frequencies and 
coupling constants with the cavity field. Here, we consider 
the symmetric collective excitation subspace V'-'l of dimen- 
sion (A^ -I- 1). We now introduce the collective operators. 



N , N 



"=1 



(3) 



which obey the following angular momentum commutation 
relations. 



[7,,7±] = +y±, [J+,J-]^2J,, 



(4) 



Furthermore, we define the Dicke basis vectors \j,m} (j - 
N/2, m = -j,-J+ 1, ■ ■ ■ ,j- 1,7), which satisfy J^\j,m} - 
jij + and J.\j,m} = m\j,m). One can conclude 

straightforwardly from Eq. (3) that the magnetic quantum 
number m equals the half difference between the numbers of 
excited atoms and the ground state atoms. In terms of the 
Dicke states, the symmetric excitation subspace, V^^^, is 



\j,m) = N,„Ji 



t=\ 



SO 



j+m 



\G), 



(5) 



According to Eq. (3), we can find 



J±\i,m) = ^crf \j,fn) 



= ^{j±m+ \)(j + m) \j, m+l). 



(6) 



It follows from Eq. (6) that, the ladder operators J± describe 
the action of pumping one more (7+) or less (/_) atom from 
the ground state \g) to the excited state \e). Accordingly, the 
magnetic quantum number m increases or decreases by one. 
Therefore, when the ladder operator 7+ acts on the collective- 
excitation state s {0 < s < N) times, there will be s atoms 
being excited, and the magnetic quantum number m will in- 
crease by s accordingly: namely \ j, m) — > |y, m + s), which is 
implied in Eq. (5). As for the ladder operator the effect 
is inverse. Therefore, the variance of the magnetic quantum 
number m represents the variance of the atomic-collective- 
excitation number 

In terms of the above collective operators, the Hamiltonian 
(2) can be rewritten as 

Hi = cja\+{(OA + w) y-+2Wfl"'fly-+y (7+7- + . (7) 

In the interaction picture defined with respect to the free 
Hamiltonian, //free = coa'a + {coa + W)Jz, the Hamiltonian 
reads 



Hf = e (ha) J, + — {J+J- + J-J+) , 



(8) 



where «„ = a'a and £(«„) - 2Wh„. The effective Hamilto- 
nian (8) is photon-number dependent. This is a special case 
of the LMG model [25] with V = 0. The LMG model can 
also be implemented using superconducting circuits [47, 48]. 
Through the relations (J+J- + J-J+)/2 - + Jy, the Hamil- 
tonian can be expressed as 



h\''> ^-W[(J,-h,f-hl-J^]. 



(9) 



As is well known, the LMG model possesses a critical point, 
at which a QPT occurs. On either side of the critical point, 
the number of excited atoms of the ground states are different; 
thus the ground states are essentially different [28, 49, 50]. 
In our system, a similar critical point also exists. To see this 
effect explicitly, we calculate the ground state for the above 
Hamiltonian in the next section. 

The last two terms of Eq. (8) describe the interaction among 
atoms induced by photons in the cavity. This interaction be- 
tween atoms is intrinsically caused by the hopping of photons 
between different atoms. And the hopping of photons induces 
a second-order indirect interaction among atoms. On account 
of this interaction, the system shows an obvious nonlinearity 
with respect to the excitation number, as shown by Eq. (9). 



C. Quantum Phase Transition Behavior of the Ground state 



where = ^J(j - my./[(2jy.(j + m)]] and |G> = 

1^1,^2, ■■■ ,gN}- 



We now analyze the discontinuous change of the ground 
state symmetry when varying the photon number For a 
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given Fock state of the field, e(«a) is a definite c numbe 
For a general photon state li/r) we replace e{ha) by its mea 
value such as si{ha)) [or e(«a)] when our studies only cor 
cern the atomic ensemble. According to Eq. (9), the eiger 
states of the system are the common eigenstates of {J~, 
{ \j, m)\m^-j,-j+l,...,j-l, j], for ; = N/2, that is. 



i/f*|;,m>|(A> = £r U'«>I'A), 



(0) I 



(K 



with eigenenergies. 



= -W [{m - riaf - nl - j {j +!)] = w„,. (11 

Clearly, the ground state is photon-number dependent, that i; 

< «„ <j-' 

- i A ^ n 

(1^ 



|G> 



I;, , na> j- 

\j, -j) or I;, j) , n„ = 
I;, -;■) , «„ > 0, 



A <0, 
A <0, 

A > 0, 
A>0, 



where [n] denotes the (half) integer nearest to n. This fac 
means that the ground state symmetry changes suddenly whe 
the photon number is varied from one domain to another. 

In the following discussions, we restrict the analysis to th 
negative detuning A < 0. As shown in Fig. 2, when th 
vahie of the photon number na is varied in the domain c 
< «n < j - 1/2, the atoms will experience different groun 
states, which implies that QPTs occur. 

There are energy-level crossings at «^ - j - n/2, (n - 1, : 
5, ■ ■ ■ ,2j-l). In the domain j-{n + 2)/2 < na < j-n/2,th 
ground state of the system is |G) - \j,j - (n + l)/2), whei 
as in the next domain j -n/2 < na < j - (« - 2)/2, the groun 
state of the system is |G) - \ j,j - {n - l)/2). If na increase 
from j - (n + 2)/2 to j - (n - 2)/2, the energy level of th 
excited state crosses the energy level of the ground state i 
na - j- n/2. At the level crossing, the excited state \ j, j - {n - 
l)/2) and the ground state \ j, j-{n + l)/2) are degenerate. On 
the right side of this critical point the original excited state 
\j,j - (« - l)/2) in the domain of j -{n + 2)/2 < na < j -n/2 
will become a new ground state for the system in the domain 
of j - n/2 < na < j — {n — 2)/2, which implies that a QPT 
occurs. In this sense, we can use the average photon number 
«a to control the occurrence of the quantum phase transition. 
At the critical point nj, - j - n/2, both |y, j - (n + l)/2) and 
\j, j - (« - l)/2) are the ground states. Moreover, at this point, 
the ground state is highly degenerate, thus the system is in a 
symmetric phase. 

In other domains, namely, when A < and n„ > y - 1 /2, 
or, A > 0, the ground state is \ j, j) or \j, - j). In these cases, all 
the atoms are fully polarized. As all the two-level atoms can 
be considered as quasispins, the system is ferromagnetic in 
this domain, and the rotational symmetry is broken. Thus the 
system is in a symmetry-broken phase. Notice that in the left 
vicinity of the critical points n^, under the condition A < 0, 
the ground state is \ j, m = [««]) and possesses one less atomic 
excitation than that in the first excited state \j, [na] + 1). It is 
clear that \j, [ng]) and \j, [na] + 1) are nearly degenerate, but 
their energies are much less than that of \ j, [na] + 2). Thus, 
there exists an energy gap that makes exciting two more atoms 
difficult, but easy for exciting one more atom. We call this 
effect "atomic blockade." 
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FIG. 2: (Color online) Diagram of the ground state of atoms con- 
sisting of A' two-level atoms controlled by the cavity photon number 
when A < 0. (a) Diagram of the energy levels versus the magnetic 
quantum number m. The upper figure in (a) shows: the ground state 
located at m = n„ = j - (n+ l)/2; The lower figure in (a) shows: the 
two degenerate ground states located at ra = j-(;i + 1)/2, j-(n- 1)/2, 
respectively, while n„ = j - n/2; (b) Diagram of the ground states 
corresponding to different average photon numbers in the cavity. 



in. DMVEN ATOMIC ENSEMBLE 

As there exists a level crossing for the photon-dressed 
atomic ensemble at na - «„, we apply a weak classical driv- 
ing to the atomic ensemble. The interaction is described by 
the Hamiltonian, 

N 

H.^^Yj (o--^'^^""' + erf e"'"''") , (13) 

where Q. is the Rabi frequency and is the driving frequency 
of the drive. The total Hamiltonian H - H\ + H2 becomes 



fj(R) ^ fjy-i + + ^ _ ^^,) + o (/_ + j^) (14) 



A') 



in a rotating frame with rotation exp[/(wj/- + toa^a)!]. In this 
driven case, the photon number a ' a still is a conserved num- 
ber. Therefore the photon number does not change in the dy- 
namical evolution even though we apply a classical driving 
field. As a result, we can treat the photon number as an in- 
dependent external parameter, which is decoupled from the 
atomic dynamics. We tune the frequency ajj to satisfy the 
resonance condition cja + W - aid =0. Then the simplified 
Hamiltonian is H^^^ = h\'^ +H' with //' = Q (J- + J+). When 
the optical field is prepared in a coherent state \a}, the Hamil- 
tonian, after this average «„ — > «„ = («o), reads 



= -W [(J, - n„f -n]- J-\ + Q.(J- + (15) 

where («„) = lap, for «„ = 1/2 + 6. Here 6 is the devia- 
tion from the degenerate (critical) point. To see if the atomic 
blockade effect occurs, we express the above averaged Hamil- 
tonian in the angular momentum basis as 



H^"^ = ^ aj,„\j,m){j,m\ 

m=-j 

+ Q,„+i (I;, m+l) {j, m\ + h.c.) , (16) 



where = Q -m + \){j + m). We can then more read- 
ily observe the transition from |j, m) to \ j, m + 2) by exciting 
two more atoms around the critical point n^. 



A. Reduced dynamics on the subspace with m = 0, 1 

When the photon number is in the vicinity of 1 /2, the 
nearly degenerate energy levels m -0,1 {\j, 0) and \j, 1)) will 
be strongly coupled with each other as a result of the driving, 
but weakly coupled with other energy levels. Then the two 
energy levels (m = 0, 1) form a relatively stable subsystem. 
Hence we can treat the transitions from the subsystem to other 
levels by a perturbative approach. In terms of the states with 
definite quantum number m, the Hamiltonian //'^^ = Hq + Hi 
can be decomposed in two parts, the nonperturbative Hamil- 
tonian, 

Ho = woUO)0;0| + a;i|y,l>0;i| 

+QVj(; + i)Ui>O;0| + h.c., (17) 

and the perturbation, 

H: = Q2|;,2)0-,l| + Qo|;,0)0-,-l| 

j 

+ ^ i^m \i, m) {j, m\ 

m=-J,mtO,l 

+ Yj n,„+il;;ni+ 1)0;ot| +h.c.. (18) 

To see clearly if the atomic blockade effect occurs, namely, 
if it is difficult to excite two more atomic excitations, we need 



{m=^Q,l-j<m<j) 



(m^O, I -j<m<j) 




(a) 



(b) 



FIG. 3: (Color online) Energy-level diagram of the m = 0, 1 subsys- 
tem of the driven atomic ensemble, (a) The two nearly degenerate en- 
ergy levels \j, 0> and \j, 1> are strongly coupled with each other, when 
the average photon number in the cavity is ii„ a 1 /2, but weakly cou- 
pled with other energy levels, {b) The effective subsystem spanned 
by \Ao) and \Ai) when using the perturbation approach. 



to find the transition amplitude for the system initially pre- 
pared in the subspace spanned by |y, 0) and \ j, I) to the doubly 
excited state \j, 2) around the critical point via - 1/2. To make 
ly, 0) and \ j, 1) nearly degenerate, we restrict Q < iia < 1. We 
note that we can also choose any other pair of nearly degener- 
ate states around the corresponding critical point which makes 
the pair nearly degenerate. We first diagonalize the nonpertur- 
bative Hamiltonian (17) as. 



//o = ^ol^o)aol + /iil^i)<^i|. 
The two eigenstates are 

\Ar)^A-'[^Aj,Q) + \jA)], r = 0,l. 



with corresponding eigenenergies. 



where A,. = -i- 1 are normalization constants with 

W6 + (-\yp 



and 



(19) 
(20) 

(21) 
(22) 
(23) 



We note that {j,m\Ar} - for m ^0,1. Therefore, \Ao}, \Ai} 
and \j, m) (m ^ 0, 1) form a complete basis of the Hilbert space 
for a given j. In this basis, Hj can be expressed as, 

H, = Q2[^iU2><^ol + /72U2)(Ai|] 

+Qo [/73l^o>0;-l| + ?74l^i)0-,-l|] 

j 

+ ^ <^m \j, m) {j, m\ 

m=-j.mi^O,l 

+ "'"+1 \J' + 1> 0'. m\ + h.c, (24) 

m=-j,mi^-l,0,l 



) 



where 



772 = - 



Ao Ai 
m = 7 T' '74 



(25) 



which satisfy |77i|- + \ri2\~ - 1 and I773I + I/74I = 1. It fol- 
lows from Eq. (24) that the transition between |/lo) and \ Ai} is 
inhibited, which is shown in Fig. 3. In order to calculate the 
correlation function g^^^ with the perturbed Hamiltonian, we 
move to the interaction picture by choosing, 

./ 

H'^ = ^0 l^o> <^ol +^1 \M 1^'' "'^ '"I ^^^^ 

m=- j,mtO, 1 

as the free Hamiltonian. In the interaction picture, the Hamil- 
tonian = //,', + H'j, where 



a,+il;,m+l)0',m|+h.c. (27) 

m=-j.iiitO 



becomes 



V,(t) = ^2 I;, 2) (7,1 <^ol + 7,2 | e''*^ '') 

+Qo (773 Mo> e-'*-'"' + 774 |ii > e-'^-'-'') {-l,j\ 
j-i 



(28) 



which is time-dependent. Here, we have defined 

A,,,',,- = OJ,,,' - /f,-, <^m,l = (^m - W/, (29) 

where 77?' 0, 1, r = 0, 1, and A„,' ,• is the energy difference 
between the diagonalized almost-degenerate energy levels la- 
beled by \Ar) (r = 0, 1) and the other energy levels labeled by 
\j,m) (m 9i 0, 1). 

B. Reduced dynamics on the subspace with m = j - l,j 

Here we consider the effect of the QPT on the higher-order 
quantum coherence around the critical point iia - j - 1/2. 
Similar to the previous section, it can be seen that the states 
\j,j - 1) and \j,j} form a relative stable subsystem. We can 
also treat the transitions from the subsystem (772 = j - l,j) to 
other energy levels by using a perturbative method. To this 
end we diagonalize the Hamiltonian in the subspace spanned 
by the two nearly degenerate energy levels \j, j - 1} and \j, j). 
It follows from Eq. (16) that, the nonperturbative Hamiltonian 
is. 



Ho 



+nj\jj} {jj -i\ + h.c. 



11 



(30) 



: \J,m) 

_j_ [m^j-Xj, -j<m<j) 



_^{m^j-Xj, -j<m<j) 



rv 



(a) 



(b) 



FIG. 4: (Color online) Energy-level diagram of the subsystem com- 
posed of m = j - I, j of the driven atomic ensemble, (a) The two 
nearly degenerate energy levels \j,j - 1> and are strongly cou- 
pled with each other, when the average photon number in the cavity 
is ria X j - 1/2, but weakly coupled with other energy levels, (b) 
The effective subsystem spanned by |/l[,> and |/l|> by the perturbative 
approach. 



with the eigenenergies, 

A',. = 2-i[-l- 277„ + 4;(l+««)]W 

-i-2-'(-l)'+y, r = 0, 1, (31) 

and eigenvectors, 

U;:) = (A^)-' [^';. \j, j-l} + \j, j)] , r = 0, 1, (32) 



where A^. - -sjl^rf + 1 (t" = 0, 1) are normalization constants 
with, 

^1. = [2Q V27]"' [(-1 + 2j - 2na)W + (-1)'-+'/] , (33) 



and 



p'= ^(1 - 2; + 277j2H^2 + 8;Q2. 



(34) 



Similar to the above Subsec. Ill A, we also note that 
{j,m\A'i.) = for 777 j - l,j. Therefore, and 777) 

{m j - I, j) form a compete basis of the Hilbert space for 
a given j. In terms of and |/lj^, the residual terms of the 
Hamiltonian (16) //; = - //(j read as 

= 0;_i [(775 I^S) + 14 1-^0) <j' > - 21 + h-c] 

j-2 

+ ^ w„, I;, 772) {j, m\ 



i«=-j 
j-i 



+ ^ Q„+i 777 + m\ + h.c.) , 



(35) 



where we used the expressions, 
\j, j) = 77^ |4) + 77^ ) , I;, - 1) = 775 |4) + 77^ , (36) 
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with coefficients defined by 



'74 



(37) 



which satisfy ^rj'^^'' + jT/^l = 1, \rj'j^\ + \rj'^^\ - 1. It follows 
from Eq. (35) that, there is no transition between and 
which is shown in Fig. 4. In order to change to the interaction 
picture, we choose the diagonalized terms, 



j-2 



H'o' = 4 |4) (4| + A\ l^-,-) (^',-| + ^ \j, m) a m\ (38) 



m=-j 



as the free Hamiltonian, and the corresponding interaction 
Hamiltonian, 



J — 



(39) 



Finally, we obtain the interaction Hamiltonian, 



J-3 



+ ^ I;, m + 1) m| e™- '"' + /i.c (40) 



m=-] 

in the interaction picture where 



(41) 



for m' 5t y - 1 , y, r = 0, 1 and AJj^, ^ is the energy difference 
between the diagonalized almost degenerate energy levels la- 
beled by (r = 0, 1) and the other energy levels labeled by 

\i,m) (m iz 

Note that Figs. 3 and 4 show transitions between three level 
systems, where some of the transitions are turned on and off. 
Indeed, it is also possible to turn on and off transitions be- 
tween three energy levels using artificial atoms made of su- 
perconducting qubits [51]. 



IV. STATISTICAL PROPERTIES OF THE ATOMIC 
EXCITATIONS 

Since /_ (7+) can decrease (increase) a single excitation, 
their roles are similar to the actions of the annihilation (cre- 
ation) operator of photons a (a ' ) for the usual bosonic system. 
Using the Holstein-PrimakofF transformation [52], the angu- 
lar momentum operators can be expressed in terms of a single 
bosonic mode. 



J+ = Va^^^^, J = ( Va^ - b'^b^ b, J- = b'^b ■ 



N 

2' 
(42) 



The angular momentum operators will become bosonic op- 
erators in the limit of large and low excitations, namely, 
(b'^'b) <*: [29]. Specifically, in this condition, one can ex- 
pand the square term VA^ - b'b on the order of {b^b)/N and 
keep to the zeroth order of b^'b/N. Then it is straightforward 
to see that J+ ^ b'^' ^^N and J ^ ^/Nb [29]. Then we can 
define a generalized second-order coherence function. 



^(2), _ {JAt)JAt + T)J-(t+T)J-{t)} 

' < A (0 J- (0) (Mt + T)JAt + r)) ' 



(43) 



for the symmetric collective excitations of the atomic ensem- 
ble, which can be regarded as a normalized correlation func- 
tion. Please note that this definition is in normal order on 
the angular momentum operators /+ and which satisfy 
that the average of any analytical function of in nor- 
mal order over the ground state \j, -j) equals zero, that is, 
0' ~j\ ■ fiJ+J ) ■ \ j^ -j) - 0. Here j is a conserved quantity. 
This property satisfies the conventional normal order defini- 
tion about the bosonic operators (0| : f{b'^b) : |0) = in the 
second coherence function. This coherence function g^~\T, t) 
is proportional to the joint probability of observing one ex- 
cited atom at time t and another one at time t + T. To study the 
generalized second-order coherence function g^^'(T, t) in the 
stationary state, below we consider it in a unitary evolution 
case (without dissipation) and also in a dissipation case but at 
a steady state. 



A. The case without dissipation 

Firstly, we consider the case where the system is free of 
dissipation. In this case, (■ ■ ■ ) in Eq. (43) for §'~'(r, t) denotes 
the average of an observable over the initial pure state. 



|iA(0))= 2c,„|;;m), 



(44) 



where YJ,„=-j\cm\^ = 1- 

We next calculate the generalized second-order coherence 
function around the point 5 - Q (i.e., iia - 1/2). Since 
U{t) = Uo(t)U,(t), where Uo(t) = exp(-///^r) and Ui(t) = 
T exp[-i Vi(T')dT'] are the free evolution and the dy- 
namics due to the interaction, respectively. We note that 
Uq{t)J+J-Uo(t) - J+J- is useful in the following calcula- 
tions. Using this result, the generalized second-order coher- 
ence function g^'' (t, 0) becomes 



{iif'mu]{T)j^j-u,{Twm 
{rmnmm\u](T)j^j^u,{T)\m) 



(45) 



where |i/''(0)) = 7-|i/'(0)). We will calculate analytically 
the generalized second-order coherence function by apply- 
ing standard perturbation theory, with Vi{t) as a perturbation. 
Let us first consider the conditions where the perturbation ap- 
proach is valid. If we tune the atom-field detuning A and the 
Rabi frequency Q. of the driving field to be suitable and place 



an appropriate number of atoms in the cavity, we can make 
the perturbation theory valid, that is, for m' = 2, -1, r = 0, 1, 
and m 9^ - 1 , 0, 1 these conditions expUcitly are 



Q2771 « A2,o, Q2?72«A2,i, 

Q0773 « A_i_o, Qo'74 A_i_i, 



(46) 



Under the above conditions, we can treat the time-evolution 
operator C//(t) perturbatively. When ria is in the vicinity of 
the critical point «JJ (for m - 0,1, n'^ - 1 /2), the energy levels 
of \Ao) and \Ai} are nearly degenerate. The energy differences 
Aij and <x>m+\,m ('« 5^ -1,0, 1) are very large compared with 
the level spacing between |/lo) and \Ax). Hence, under this con- 
straint, the above conditions (46) can be satisfied by varying 
the Rabi frequency Q.. Since the state \ j, 0) is the ground state 
when < < 1/2, then \ j, 1) is the state by exciting one 
more atom. Similarly, \ j, 2) has two more excitations than the 
ground state, and has a much higher energy than that of \j, 1). 
However, \j, 1) is the ground state when 1/2 < < 1, yet 
\i, 0) is an excited state which has one less atomic excitation 
than the ground state \ j, 1). To see if two excitations are sup- 
pressed, we choose cq - c\ - \ I V2 and c,„ = (for m 5^ 0, 1) 
in the initial state. 



1 



|iA(0)> = col7,0> + cil;, 1) = — (|j,0> + \j, 1». 

V2 



(47) 



When the average photon number is in the vicinity of njj — 
1/2, the states |y, 0) and |y, 1) are nearly degenerate. No- 
tice that here the average photon number is in the domain of 
< n„ < 1 and around - 1/2, that is, -1/2 < 5 < 1/2, 
and \6\ is very small. Using first-order perturbation theory, 
the generalized second-order coherence function in Eq. (45) 
is approximately 



X 



where 



X = Xi -H X2 -I- X3 -I- X4 -I- X5, 



(48) 



(49) 



The parameters Xf, (^1 - 1, 2, ■ ■ ■ , 5) and jf, (£2 - 1,2,3) 
have complicated expressions, which are presented in the ap- 
pendix. The generalized second-order coherence function 
given by Eq. (48) is illustrated in Fig. 5. It is shown that, as 
increases, the value of g^^^r, 0) approaches unity with some 
oscillations. Physically, Eq. (48) describes the joint probabil- 
ity of observing one excited atom at instant f = and another 
after a time interval r. In Sec. V, we use Eq. (48) to analyze 
the atomic blockade effect. 



B. the case with dissipation 

In this subsection, we consider the system surrounded by 
a thermal reservoir at zero temperature. When the system is 



1 
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FIG. 5: (Color online) Second-order correlation function g'^'(T, 0) 
versus the time interval r for A' = 2 (blue curve), N = 4 (red curve), 
N = 6 (cyan curve), and A' = 10 (magenta curve), respectively, in 
the case without dissipation. Recall that g'^^(T, 0) is proportional to 
the joint probability of observing one excited atom at time t = and 
another one at time t. The first-order approximate results are shown 
using dashed curves and the exact numerical results are shown using 
solid curves. The dashed curves overlap with the solid curves very 
well. Other parameters are go = 100, go/A = -0.1, QN = 0.1|W|, 
S = -0.02. 



prepared in a state with density operator ps, the generalized 
second-order coherence function is written explicitly as 



Tr[J^MT)J4T)J-pAt)] 
Tr[7+7-p,(f)]Tr[7+y-p,(f + t)] ' 



(50) 



According to Eq. (50), we need to calculate the time- 
dependent density operator p.s(f) of the atoms. In the regime 
of weak coupling of the driving field [53], which demands the 
driving field to only perturbatively change the energy levels, 
and assuming the atomic ensemble to be in a common reser- 
voir, then the master equation is approximately 



dps (t) 
dt 



J-Ps it)J+-l{J+J-,Ps{t)} 



(51) 

where y is the collective decay rate of the atomic ensemble. 
Since the photon number is a conserved number, and the fre- 
quency of photons is in large detuning, it does not influence 
the dynamical evolution of the atoms. Then the influence of 
the decay of the photons is negligibly small to the atoms. We 
resort to numerical calculations to show the results about the 
steady state by plotting g'^'(0, f — > 00) versus 3 in Fig. 6 (a) 
and g^^'(T, t — > 00) versus r in Figs. 6 (b)-(d). By comparing 
them with the results in Fig. 5, we will discuss them in the 
next section. 
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V. DOUBLE EXCITATION EFFECTS I: THE ATOMIC 
BLOCKADE EFFECT 

In this and the next section, we discuss some physical ef- 
fects due to the double collective excitation, according to their 
quantum statistics characterized by the generalized second- 
order coherence function g^'^r, t) introduced in the last sec- 
tion. We have calculated the generalized second-order coher- 
ence function in the above section both in the dissipation-free 
case and also the case with dissipation. In this section, we 
discuss the results in both cases according to the above cal- 
culations. We illustrate the analytical results (48) and com- 
pare them with the numerical results by plotting in Fig. 5 the 
generalized second-order coherence function g'~'(T, 0) versus 
the time interval r around 6c - 0, without dissipation. The 
generalized second-order coherence function is plotted for 

= 2, 4, 6, 10 atoms, respectively. It is clear from Fig. 5 that, 
close to the critical point 6c - 0, our analytical approximate 
results (48) (dashed line) agree very well with the numerical 
result (43) (solid line). Obviously, g^^\T, 0) < 1 at any time 
interval t. This shows that the atomic collective symmetric ex- 
citations obey sub-Poissonian statistics. It can also be found 
that as increases, g*^^(r, 0) < 1 oscillates slower and slower 
and approaches unity, especially for g^~*(0, 0). That is because 
the generalized second-order coherence function at r = is 



g<-*(0,0) = 1- 



N- + 2N 



(52) 



Hence (0, 0) increases as increases. In the thermody- 
namic limit — > oo. 



(53) 



This shows that when A^ is smaller, the effect of suppressing 
the doubly-excited state is enhanced. 

Figure 6 shows the results for the dissipative case. Fig- 
ure 6(a) shows g'^'(0, oo) versus the average photon number 
6 in steady state for A^ = 2,3,5, 10 atoms, respectively. As 
shown in this figure, the value of g^~\0, oo) increases as A^ in- 
creases for a larger average photon number For a definite A^ 
and a small value of 6, g^^\0, oo) increases as 6 increases. At 
some intermediate time there is a peak in g^^\0, oo) followed 
by a steady decrease, asymptotically approaching a constant 
value for large 6. The smallest value of g^~\0, oo) occurs at 
6 = -0.5. This phenomenon is also prominent in Figs. 6(b)- 
(d). Figures 6(b)-(d) show ^*^'(t, oo) versus t for N = 5 and 
6 = -0.5, and j - 1.1, respectively. The antibunching effect 
of collective excitations of an atomic ensemble is observed 
since the envelop of g*^'(T, oo) shows g^^\Q, oo) < g*-'(r, oo) 
with some increasingly rapid oscillations as 6 increases in 
Figs. 6(b)-(d). Additionally we note Hiat g'^Xr, oo) approaches 
one, as expected, after some oscillations. This indicates the 
probability of two collective excitations of the atomic ensem- 
ble at the same time (r - 0) is smaller than that at a different 
time (r 0). Therefore, the resonant excitations from the 
ground state to the doubly excited state are suppressed. This 
is a clear signature of the atomic blockade. Compared with 
Fig. 5, this result is better and closer to physical reality. As 




0.6 .' 



0.5,' 

0.4- 




(b) 5=-0.5 




FIG. 6: (Color online) Numerical results for the generalized second- 
order coherence function g*^'(0, oo) and g'^*(T, oo) with dissipation in 
the steady state, (a) g*^'(0, oo) versus 6 for N = 2 (blue thick solid 
curve), N = 3 (red dashed curve), N = 5 (cyan dashed-dotted curve) 
and N = 10 (black thin solid curve), respectively; (b)-(d) g*^'(T, oo) 
versus t with A' = 5 for <5 = -0.5, 6 = 0, and 6 = J - 1 . 1 , respectively. 
Other common parameters are y = 1, go = 100, go/A = -0.1, and 
flN = OA\W\. 



shown in Fig. 5, the generalized second-order coherence func- 
tion only oscillates with time interval t and does not approach 
1 as we expect when r ^ oo. In Ref. [54], the photon anti- 
bunching eff'ect is also obtained in only two interacting atoms. 
However, the antibunching effect we obtain is about atomic 
collective excitations, and the photon number is a conserved 
number. In this sense, we do not need to consider the photon 
correlation. 

To conclude this section, we give some remarks about the 
atomic blockade. For applications in quantum information, 
the atomic blockade provides a novel approach to physical im- 



plementation of scalable quantum logic gates such as imple- 
menting a CNOT gate between two atoms [11-13] and some 
kinds of quantum protocols [55-58]. Furthermore, as double 
excitation are inhibited in the Rydberg blockade mechanism, 
it also supplies a fascinating approach to store quantum infor- 
mation [11, 12]. However, the dipole-dipole interaction de- 
pends on the distance between Rydberg atoms. To achieve a 
stronger interaction, it requires the atoms to be closer in space 
or to be excited to higher Rydberg states, in which the prin- 
cipal quantum number is very large, but this will not be con- 
venient to control the atoms individually [11-13]. Such as 
the Rydberg levels n = 79 and 90, the corresponding block- 
ade shift is 27r X 3 and 2n x 9.5 MHz at an interatom dis- 
tance 10.2 jum, respectively. To achieve a larger energy-level 
shift due to the Rydberg blockade, the distance needs to be 
decreased, and thus the coherent manipulation of individual 
atoms is difficult. It is this consideration that motivates us 
to find a new mechanism inducing a stronger interatom cou- 
pling, valid for long distances and controllable to improve the 
dipole-dipole interaction. We note that in Ref. [59], the cou- 
pling strength between atom and photons can reach 2;7!- x 120 
MHz in a high-finesse cavity, which leads us to anticipate that 
the strong atom-photon coupling will induce a stronger inter- 
atom interaction among atoms. In addition, this interaction 
can be feasibly controlled by the volume of high-finesse mi- 
crocavities. This fact means that to achieve a strong interatom 
interaction among atoms will not take stringent requirements 
on manipulating atoms individually. Therefore, from the point 
of view of the controllability and strength of the interaction, 
the photon-induced interaction among atoms in our system is 
better than the dipole-dipole interaction inducing the Rydberg 
blockade. 



VI. DOUBLE EXCITATION EFFECTS II: SENSITIVITY 
OF THE QUANTUM PHASE TRANSITION 




3.5 3.6 3.7 3.8 3.9 4 4.1 4.2 4.3 4.4 4.5 

5 




8.5 8.6 8.7 8.8 8.9 9 9.1 9.2 9.3 9.4 9.5 

5 



FIG. 7: (Color online) Numerical results for the generalized second- 
order coherence function g*^*(T, 0). Here: f = 0, r = 3, go = 1000 
(red thick solid curve), go = 3000 (black thin solid curve), (a) N = 
10, Sc = 4, (b) N = 20, 6c = 9; other parameters are the same as in 
Fig. 5. 



As the system possesses QPTs, we now analyze how to con- 
trol the QPT by photons in the cavity. To show the effect of 
the QPT on (r, 0) more clearly, we consider the g*-* (r, 0) 
around the critical point = j - 1/2 at a fixed time interval 
T. Then, according to Eq. (44) we choose cj-i - cj - 1/ V2 
and c,„ = (m 5t j - 1, j) in the initial state, namely 



\fm-^\jj-i} + ^\jjy 

V2 V2 



(54) 



With the relations U (t) = U^^(T)U'j(T) for U'^(t) = 
expiiH^'r), Vjir) = T exp[-/ j^'" y;(T')c/T'], it follows from 
Eq. (43) that 

UI\t)J^J^U'q(t) 
= /a(f)(|;;;-l)<;,;|-|;,;)0;;-l|) 

+Af)(l;,;-i>0-,;1 + l;,7)0;7-i|) 
+ymi,mj\-\j,j-i){jj-^) 

+J+J-. (55) 

The explicit expressions of the coefficients a (f), y6(0 and y (f) 
are given in the appendix. 



Next, we use the perturbation approach to calculate the gen- 
eralized second-order coherence function under the following 
conditions for m' = j - 2, r - 0,1, and m j - 2, j - I, j. 



n 



'm+l 



1)1+ 1 ,m ■ 



Up to first order in Vj (t), we obtain 



(37-i)(i:Li4)' 



(56) 



(57) 



where the parameters jc^^ {{i - 1 , 2, ■ ■ ■ , 7) and y^^ ((2 - 
1, 2, 3, 4) have very long expressions, so we give these in the 
appendix. 

We also numerically calculate the generalized second-order 
coherence function in Eq. (43), and then plot g^^\T, 0) ver- 
sus 5 in Fig. 7. As Fig. 7 indicates, the statistical coher- 
ence of atomic excitations is very sensitive to the critical point 
rfg— j - 1 12. The probability of double atomic excitation is 
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above the dotted straight line in the left-hand side of the crit- 
ical point, while it is below this curve in the right-hand side 
of the critical point. Furthermore the envelope exponentially 
decays. When the average photon number is in the domain 
of y - 1 < «fl < j - 1/2 (j - 3/2 < 6 < j - I), the energy level 
of \j, j) is higher than \j, j - 1) but lower than \j, j - 2); while 
in the domain of [y - 1/2, j], the energy level of \j, j) is lower 
than both \J, j - 1} and |j, j - 2), and the order of the energy 
levels is ojj < coj^x < Wj_2 < < We also note 

that, as the coupling strength increases, '(t, 0) oscillates 
faster with respect to 5. In addition, as increases, the value 
of g^~*(r, 0) increases. 

Above, we gave a qualitative analysis of the generalized 
second-order coherence function based on perturbation the- 
ory. According to our calculations, there is a large dis- 
crepancy between the theoretical analysis and the exact nu- 
merical result. The reason may be as follows. As seen 
in the definition of the generalized second-order coherence 
function, [i.e., Eq. (43)], this is determined by two corre- 
lation functions, that is, ((/r'(0)|f/f'"(T)7+7-t/;(T)|tA'(0)) and 

<iA(0)|f/f"(r)y+y_f/;(T)|i/r(0)>. As far as the latter is con- 
cerned, we calculate it in the interaction picture. Here, we 
approximate the time-dependent wave function [/y(T)|i/r(0)) to 
first-order by perturbation theory. Since the operator J+J- 
gives two large and markedly different eigenvalues to the com- 
ponents \i, j) and \j, the originally small deviation in the 
approximate wave function with respect to the exact one will 
be enlarged. 

However, when we come to the case with m = and 1, 
the situation turns out to be totally different. First of all, let 
us turn to the Hamiltonian H - Hq + Hi given in Eqs. (17) 
and (8). In the large-detuning regime, it only induces a Rabi 
oscillation between the two nearly degenerate states \ j, 0) and 
\j, 1), while leaving the populations in the other states almost 
unchanged. On account of the conservation of the total prob- 
ability and the same eigenvalues of the operator J+J- on the 
two relevant states, in the system which is initially in an equal 
superposition of |y, 0) and \ j, 1), the approximate correlation 
function (i/'(0)|t/^(T)7+7_f//(T)|i/'(0)) is expected to be quite 
close to the exact one. This situation will not take place for 
the case with m - j and j — I, since the relevant eigenvalues 
of the operator J+J^ are remarkably different from each other. 
A similar analysis can be applied to the numerator in the gen- 
eralized second-order coherence function. Consequently, the 
generalized second-order coherence function obtained from 
the perturbation theory will coincide with the exact one for 
the case with m - and 1, while there is an obvious differ- 
ence between the results from these two methods for the case 
with m - j and j - 1 . Therefore, we only give the numerical 
results in Fig. 7. 



VII. CONCLUSION AND REMARKS 

In this paper, we study the statistical properties of atomic 
excitations for two cases: with dissipation and without dis- 
sipation. We find that this statistical property can be con- 
trolled by the average photon number in the cavity. In ad- 



dition, the photon-induced second-order interaction between 
atoms is valid in the long range and can be strengthened by 
a high-finesse microcavity with a very small effective mode 
volume. Furthermore, we find that the double atomic excita- 
tion will be suppressed when the average photon number in 
the cavity is in the vicinity of some special points (degenerate 
points). We have also studied the critical behavior of this sta- 
tistical property of atomic excitations around the critical point 
at which the QPT occurs. 

To characterize the statistical property of atomic excita- 
tions, we define a generalized second-order coherence func- 
tion similar to the second-order coherence function for pho- 
tons. Furthermore, in the limit of — > oo and low excitations, 
it becomes the conventional one. We have demonstrated the 
antibunching effect for atomic excitations near the degener- 
ate points and the characteristic of sub-Poissonian statistics, 
which implies the existence of the atomic excitation blockade. 
Since this system possesses several critical points, we also 
study the critical behavior of the generalized second-order 
coherence function of atomic excitations around the critical 
points. Our results show the sensitivity of the system dynam- 
ics with the average photon number in the cavity. 
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Appendix: EXPLICIT EXPRESSIONS FOR THE 
PARAMETERS OF g<-' 



In this appendix, we present the expressions for the param- 
eters used in Eqs. (48) and (57), respectively. 

For m = 0, 1, the parameters of g^^^T, 0) given by Eq. (48) 
are 



Xl 

X2 



I'oO' - 1)^0' + 2)-|ci|- 1771/7302,0 + mmOijf , 



mm 



\0\o-OUa\ 



X3 = fU+lf 



mm - mm 

voC() (774771(9* j j - 773772(9* j o) 



+ Cl 



772 ^73 - ?7i'74 

X4 = 7'(; + 2)(/ - 1) |voCi (7730-1,0 + 77^C_i,i) + Co 

-^5 = via - 1)(/ - 4)(; + 3)4^ |1 - e'"--f , 

^-2,-1 



(A.l) 
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and 



where 



yi = O'+iU 

y2 = Q}{j-lf{j + 2f 

X |?7i(co'73 + ci77i)(92,o + micom + cimWxi] , 

y3 = vl(j-i)(j + 2) 

X |773 (co?73 + ci77i) (9-1,0 + J]4 (coT]4 + ci772)0-i,i| , 

(A.2) 



where 



(A.3) and 



For m - j-l,j, the parameters of g^^\T, 0) given by Eq. (57) 
are listed as follows: 

= 2;|flo(T)P, 

4 = 2(2;-l)|fli(T)|-, 

x^ = 3(2;-2)|fl2(T)P, 

4 = 4(2;-3)|fl3(T)|-, 

Xg = -2a (t) Im [ao (t) a\ (t)] , 

x^ = 2;8(T)Re[flo(T)fl;(T)], 

4 = rW[|flo(T)|2-|fli(T)P], 



and 



= i, /,^2j-h 
yl = 3(2;-2)|c2|% y^^ySW, 



(A.4) 



(A.5) 



with 



a2(T) = V2; - 1 



1 + ^/..(r) 



2;-l 



fl3(T) = nV3(2;-l)(2;-2) 



(1 -e'^^j-'-'-^O 



C2(t) 



fhiir) 

V2(2;-l)' 



«(0 = ' sin((7f), 

y6(f) = qoq'^{tOj-i - cJj) [cos{qt) - 1] , 

r(f) = 2qoq-^n^j2j[cosiqt)-l], 



- on Fr-<- n V2Q(2;-1) 
qQ = -2Q V2;(; - 1), / = — 

I 1 - e' 

q = ^(coj-i - Ujf + 8;Q2, o;„,„ s — — 



(A.6) 



(A.7) 



■.(A.8) 



Here, 
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